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In view of the application to supernova simulations, we calculate neutrino reaction rates with 
nucleons via the neutral and charged currents in the supernova core in the relativistic random phase 
approximation (RPA) and study their effects on the opacity of the supernova core. The formulation 
is based on the Lagrangian employed in the calculation of nuclear equation of state (EOS) in the 
relativistic mean field theory (RMF). The nonlinear meson terms are treated appropriately so that 
the consistency of the density correlation derived in RPA with the thermodynamic derivative ob- 
tained from EOS by RMF is satisfied in the static and long wave length limit. We employ pion and 
rho meson exchange interactions together with the phenomenological Landau - Migdal parameters 
for the isospin-dependent nuclear interactions. We find that both the charged and neutral current 
reaction rates are suppressed from the standard Bruenn's approximate formula considerably in the 
high density regime (pb ^ 10 14 g/cm 3 with pi, the baryonic density). In the low density regime 
(Pb J; 10 14 g/cm 3 ), on the other hand, the vector current contribution to the neutrino - nucleon 
scattering rate is enhanced in the vicinity of the boundary of the liquid - gas phase transition, 
while the other contributions are moderately suppressed there also. In the high temperature regime 
(T > 40MeV with T the temperature) or in the regime where electrons have a large chemical poten- 
tial, the latter of which is important only for the electron capture process and its inverse process, 
the recoil of nucleons cannot be neglected and further reduces the reaction rates with respect to the 
standard approximate formula which discards any energy transfer in the processes. These issues 
could have a great impact on the neutrino heating mechanism of collapse-driven supernovae. 



I. INTRODUCTION 



Calculations of neutrino - nucleon reaction rates in a hot (T < 50MeV) and dense (10 13 g/cm 3 < p b < 10 15 g/cm 3 ) 
core of a collapse-driven supernova are complicated problems (see, e.g., [jl]-[| for the standard rates). When the density 
reaches p b ~ 10 13 g/cm 3 , the average separation of nucleons d ~ 6x 10~ cm (p b / 10 13 g/cm 3 ) -1 / 3 becomes of the same 
order as the typical neutrino compton wave length A„ ~ 6x 10~ 13 cm (E v / 30MeV) _1 . This means that neutrinos are 
interacting simultaneously with multiple nucleons for this density or higher. If nucleons are distributed uniformly in 
space and time, which is unlikely, the outgoing waves from multiple targets interfere with one another and the cross 
sections remain the product of the number of targets and the cross section of a single scattering. However, in reality, 
the distributions of nucleons are fluctuating due to mutual interactions and the reaction rates will be modified from 
those obtained with this simple formula. Furthermore, the typical energy E v ~ 30MeV (T/lOMeV) of a neutrino, 
which is approximately an inverse of the duration time of interaction between a neutrino and a nucleon, is of the 
same order as the nucleon - nucleon scattering rate T ~ (av)n b ~ 30MeV for the density p b ~ 3 x 10 13 g/cm 3 and 
temperature T ~ lOMeV, where a, v and rib are the scattering cross section, the nucleon velocity and the nucleon 
number density. This implies that while a neutrino is interacting with a nucleon, the target nucleon is scattered off 
by another nucleon. Hence both the spatial and temporal correlations of nucleons are important. 

It is, however, very difficult to calculate these correlations of nucleons accurately, since they are induced by the 
strong nuclear interactions and cannot be treated with a simple perturbation. This issue was studied by several 
authors from early on. Sawyer [|J, for example, calculated the density and spin-density fluctuations of nuclear matter 
in the static and long wave length limit from an equation of state, while Iwamoto & Pethick [pi investigated them 
from the Fermi liquid theory. These results have not been incorporated extensively in simulations of supernovae or 
proto neutron stars (see, however, ||) and the simpler formula mentioned above of the reaction rates has been used 
even in the most sophisticated computations. |7|-|lC|l 

With increasing recognition that the neutrino transport is one of the key factors for a successful supernova explosion 
jnlO], this issue has recently attracted attentions of supernova researchers. In the neutrino heating scenario of 
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supernova explosion |]J|,[l4|], which is supposed to be the most promising at present, the shock wave stagnated in the 
iron core after the core bounce is expected to be revigorized by neutrinos diffusing out of the proto neutron star. It has 
been shown that it is important in this mechanism to increase the neutrino heating rate via the dominant processes 
of electron neutrino absorption on neutrons and electron antineutrino absorption on protons behind the shock wave. 
This heating rate is approximately given by: 
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(1.1) 



Here Y n = n n /rib and Y p = n p /rib are the number fractions of free neutrons and protons, respectiv ely; the normaliza- 
tion with the baryonic number density rib indicates that the rate per baryon is calculated in Eq. ([O]); L„,52 denotes 
the neutrino luminosity in units of 10 52 erg/s, E v i 5 the neutrino energy normalized by 15MeV, the radius in 10 7 cm; 
(fi) is the mean value of the cosine of the angle of neutrino propagation relative to the radial direction. It is clear from 
this equation that the higher neutrino energy and/or the greater neutrino luminosity increase the neutrino heating 
rate. One obvious way to achieve that is to decrease the neutrino opacities in the supernova core. Since the dominant 
source of opacity for neutrinos is the neutrino - nucleon reactions, the modifications of their rates could have a great 
impact on the supernova mechanism ]l5|-p0|]. 

As stated above, the neutrino - nucleon reaction rates used in the supernova simulations thus far were mostly 
evaluated by multiplying the rate for a single target nucleon with the nucleon number density (see, e.g., P]2l|] 
and references therein), thus ignoring the correlations due to ambient nucleons. Recently, however, some authors 
|L8|-p(lf22 23 1 calculated the density as well as the spin-density correlations of nucleons due to nuclear forces based on 



RPA. As shown later, this method evaluates the nonuniform nucleon distributions in space and time using the mean 
field approximation. On the other hand, Raffelt and his company [^5|Jl^j2^ -p6| insisted that the effect from collisions 
of two nucleons on the nucleon spin-density fluctuations cannot be neglected, which is not taken into account in RPA 
(see also |p7| ). 

One of the authors (H. T.) recently published the nuclear EOS based on RMF and the Thomas - Fermi approx- 
imation for finite nuclei covering the wide range of density, temperature and electron fraction of relevance to the 
supernova simulations p8| , p9| . In this paper, we calculate the nuclear correlations in the relativistic RPA based on 
the Lagrangian used in the calculation of EOS in RMF US 

It is shown in the next section that this guarantees 
the thcrmodynamical consistency of the neutral vector current part of reaction rates with EOS in the static and long 
wave length limit if the nonlinear meson terms are appropriately treated. The neutral axial vector current and the 
charged current contributions to reaction rates are calculated on the same basis but with additional introduction of the 
phenomenological Landau - Migdal parameters for the isovector channel 1 32 34|,|35|,^7j , since they do not contribute to 
RMF and are neglected in the theory. The possible collisional effects which are supposed to be important particularly 
in the low density regime and might have some roles in forming the neutrino spectra will be discussed elsewhere [ |36| , 
since we consider the consistency of the reaction rates with EOS which we have currently at hand is more important. 

This paper is organized as follows. In the next section we formulate the neutrino reaction rates with nucleons by 
the so-called dynamical structure functions and show that RPA is consistent with EOS obtained by the mean field 
theory. Then we represent some modifications of the reaction rates due to the RPA correlations, using the Bruenn's 
standard aprroximation formula as a reference. We summarize the paper with some discussions of implications for 
the supernova simulations in the last section. 



II. NEUTRINO - NUCLEON REACTION RATES 



First we express the reaction rates in a general form. Since we are interested in low energy reactions (E <C M w , 
the mass of weak boson), the weak interaction is well approximated by the interaction Lagrangian density: 

J C I (x) = ^l„(x)J^(x) , (2.1) 

where Gf is the Fermi coupling constant, l/j,(x) is the lepton weak current given by 

l^(x)=^ l (x)r(l- l5 )ij l (x) , (2-2) 

and Jjij(x) is the nucleon counter part, 

Jn( x ) — ^n{ x ) l^{hv - fiAjs) iPn(x) . (2.3) 
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hy and Ka are the vector and axial vector coupling constants, respectively, and are taken for the charged current 
as hy — gy = 1 and Ha — 9a — 1-23. For the neutral current they are h v = —1/2, h\ = — 1/2 qa and hy = 
1/2 — 2 sin 2 9\y, h p A = 1/2 g a for neutron and proton, respectively. 9w is the Weinberg angle. 

Following the standard procedure (see, e.g., pTj), the reaction rates are obtained by making a square of each matrix 
element evaluated up to the lowest order of the Fermi coupling constant, taking the thermal ensemble average for the 
initial state and summing over the final states: 

R(q in ,q out ) = ^ K a(3 (q in ,q out ) Sff(k) . (2.4) 

q m and q out are the four momenta of the incident and outgoing leptons, respectively, k = q m — q out is the four 
momentum transferred from lepton to nucleon. In the above equation, the tensor K a p comes from the kinematics of 
leptons and given by 

IS { „i n „ ou t\ o ( „out in i „out in „out in „ ■ • _ #7 „out in\ ( c\ r\ 

J<a0{q ,q ) = &{q a qp + q (j q a - q ■ q g a fi + « £ Q/3 q s <2 7 J • (2-5) 
Here g a /s is the metric tensor with the signature of [H ] , and e a P s ^ is the antisymmetric tensor with e 0123 = 1 . 
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All information of nucleons is contained in the so called dynamical structure function S^f defined as: 



S7{k) = ld*xe l «* (J%(x)J*(0)) , (2.6) 

where ( • • • ) stands for the thermal ensemble average of the argument. 

The structure function can be generally decomposed as follows due to the isotropy of the system: 

S% P (k) = R^k)^^ 3 + R 2 (k)(u a u' 3 - g af3 ) + R^{k)k a k p 

+ R A {k){k a u + u a k f) ) + i R 5 (k) e af3S ~< u 5Ul , (2.7) 



where u a is a four velocity of the system. Putting Eqs. (2.5) and ( p.7| ) into Eq. ( |2.4[ ), we get 

R(E m , E out , cos 9) = AG 2 F E m E out [R^k) (1 + cos0) + R 2 (k) (3 - cosf?) 

- 2(^ m +^ otlt )i? 5 (fc)(l-cos6»)] , (2.8) 

where E m and E out are the energies of the incident and outgoing leptons, and 9 is the angle between the incident and 
outgoing three momenta. There are three terms with different angular dependences. In general the third contribution 
is much smaller than the other two terms and ignored in the following discussions. The meanings of them become 
clearer if we take the non-relativistic limit for the nucleon kinematics. In this limit Ri(k) and R2(k) are reduced to 

i?i(fc)« h 2 y Jd A xe lkx { PN (x)p N (Q)) , (2.9) 
fla(fc)* ^ Jd 4 xe lkx (s l N (x)s* N (0)} . (2.10) 

Thus i?i (k) comes mainly from the vector currrent part of the nucleon weak current and is nothing but a density 
correlation function of nucleons. R 2 {k) originates from the axial vector current of nucleon and represents a spin- 
density correlation function. pn(x) is the nucleon density and s l N (x) is the spin density and the spin components i 
are summed up. 

The calculation of the reaction rates is thus reduced to the evaluation of these correlation functions. It is, however, 
easier in the field theory to treat the time ordered product instead of the ordinary one since we can apply the 
perturbation theory more easily for the former (see, e.g., |38| ). Hence we define the time ordered product of the weak 
current U^f (k) corresponding to the dynamical structure function as: 

aP(U\ — / A^ Jkx j rp jot ( ™,\ jP 



iIL a N ^k) = I dhe«* (TJ«(x)j£(0)> , (2.11) 

where the symbol T stands for the chronological ordering of operators. It is also convenient to consider the corre- 
sponding retarded Green function defined as 

m^(k) = /Ae fa e(t)([j;(4J>)]> , (2.i2) 
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where [A, B] denotes the commutator of A and B and <d(t) is the Heaviside function. In fact, the dynamical structure 
function is related with the imaginary parts of those Green functions via a simple formula: 



a0 



l + eHSCfa-A/i) ImU N 



!_ e -/3(fco-A M ) ImU N 



Ra/3 



(2.13) 



Here A/z = fi ou t — fJ-in is the difference of the che mical pote ntial s between the outgoing and incident nucleons. Thus 
what we have to do is to somehow calculate Eqs. ( 2.11) or ([2.12]). 



A. Thermodynamic Consistency 

In this section, we discuss the thermodynamic consistency of the approximation for reaction rates with that for EOS. 
By consistency we mean the following: The reaction rates are nothing but correlation functions as stated above. For 
example, the density correlation function reduces in the static and long wave length to ((-/V— (N)) ■ (N — (N))) / V 2 — 
((N 2 ) — (N) 2 ) I V 2 . Here iV and V are the baryonic number and the volume of the system, respectively. This thermal 
ensemble average is related to the thermodynamic derivative of the number density with respect to the chemical 
potential as {dN/dfi) T = {(N - (N)) ■ (N - (N)) with = 1/T, which is obtained from EOS M. This sort of 
relation should be satisfied also in approximate formulations, since it guarantees the correct behavior of reaction rates 
in this limit. In the following we show that RPA is consistent in this sense with the mean field theory fl8| . The 
argument is mainly indebted to the papers ]40|j4l[ |. For simplicity we consider the non-relativistic density correlation 
in the imaginary time formalism. The extension to the relativistic case [ fl2"f and/or to the real time formalism is 
possible p3[ |. 

EOS can be calculated once the number density is obtained for a given temperature and chemical potential. The 
number density in turn is obtained from the single particle Green function as nt(xi) = — i G(l, 1 + ). Here the Green 
function is defined as usual, iG(l,2) = (T (f)(xx)4>^ (xz)) , and 1 + in the argument denotes the limit t% — > t% +0. Hence 
the approximation for EOS can be regarded as the approximation of the Green function, or that of the self energy E 
since the Green function is determined by the Dyson equation G _1 = Gq 1 — E, where Go is the free propagator. 

The response of the system to the external disturbance is studied by adding the extra interaction like, 

£ext = Jo % ^ dx\dt\dx?,dti 4>^(x2,t2) U(x2, 35i, t\) <f>(xx,ti), where only the time integration region is shown explic- 
itly. Then the response function is given by 



6G(1,1';U) 



5U(2',2) 



[G 2 (1,1',2,2')-G(1,1')G(2,2')] = ii(l,l',2,2') , (2.14) 



u=o 



where G2(l, 1', 2, 2') is the two particle Green function, and G(l, 2; U) is the single particle Green function under the 
external disturbance and is written in the interaction representation as 

jG(12;tr) . gj^M (2 . 15) 



with S = exp(— iCext )- Eq. ( 2.14 ) is derived directly by taking a derivative of Eq. ( [2.15 ). The retarded Green function 



defined by Eq. (2.12) is obtained from the response function L(l,l + ,2,2 + ) by the analytic continuation The 



equation for the response function L(l, 1', 2, 2') is obtained by taking a functional derivative of the Dyson equation as 
L(l,l',2,2') = L (l,l',2,2') - [dx 3 dx 3 ,dx 4 dx A < L (l, 1', 3, 3') • E(3',3,4',4) • L(4,4',2,2') (2.16) 

with 

ii„(l,l',2,2') = G(1,2')G(2,1') (2.17) 

iH < 1 ' *•*•*> - Ism ■ <2 - 18 ' 

It was shown |^] that if the self-energy E is derived from a potential as E(l', 1) = 63>/6G(l, 1'), or the integrability 
condition <5E(1', 1)/SG(2, 2') = <5E(2', 2)/SG(l, 1') is satisfied, then the approximation is conservative, that is, the 
approximate number density as well as the response function so obtained satisfies the conservation equation just like 
the exact one. This result in turn dictates the behavior of the response function in the static and long wave length 
limit as follows. The response function is in general given by 
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L(fi,fe) 



1 

2^ 



dk r 



L>(k ,k) - L<(k Q ,k) 1 
Q — ko 27r 



dko 



(l„ e -/3fco) i >( fco;fc ) 

17 — fco 



(2.19) 



Here L(£l, k) is a Fourier transform of L(l, 1 + , 2, 2 + ) with respect to xi — x 2 and is an imaginary energy. L > (fco, fe) 
and L < (fco J fc) are analytic functions with respect to ko and are Fourier transformed from L > [x\, x^) — ((nb(xi) — 
(n h (xi))) ■ (n b (x 2 ) - (nb{x 2 )))) and L < (xi,x 2 ) = ((n b (x 2 ) - (n b (x 2 ))) ■ (n b (xi) - (n b (xi)})}, respectivel y. T hus 
i > (xi, x 2 ) is nothing but a density correlation function we are seeking for. In fact, we use this equation ( [2.19 ) to 
derive the density correlation function L(Q,k). If the approximate response function satisfies the conservation law, 
L > (fco I 0) = 27T(5(fco) i > (0, 0) in the long wave length limit, |fc| —> 0. This then results in the equation 



lim L(0,k) 

|h|-0 



/3L>(0,0) 



/3V((n b -(n b }) ■ (n b -(n b ))) 



(2.20) 



On the other hand, the thermodynamical derivative, (dn b /dp) T 
equation as 



-i dG(l, l + )/dp, is evaluated from the Dyson 



«9G(1,1') 
dp 



dx 2 dx 2 , G(l,2') • d ° ' (2/,2) • G(2,l') 



dp 
5(2' - 2) 



<9E(2',2) 



G(2,l') 



= - Jdx 2 dx 2 , G(l,2') 

= - Jdx 2 , G(l,2') • G(2',l') 

+ Jdx 2 dx 2 ,dx 3 dx y G(l,2') • i 3(2', 2, 3', 3) • G(2, 1' 

Solving this equation using Eqs. ( |2.16| ), ( [2.20 ), we get 

<9G(1,1+) 



«9G(3,3') 
dp 



ON 
dp 



iV '■ 



dp 



V lim L(0,k) 

fc|->0 



P((N-(N))(N-(N))) 



(2.21) 



(2.22) 



This is the equation which represents the thermodynamical consistency between the response function and EOS. 
So far the only assumption is that we have used the same approximate self energy for EOS and the response function. 
Hence the above argument can be applied to any conserved current and conserving approximation. In the next section, 
we consider a specific approximation of the self energy which leads to the mean field theory and RPA. 



B. RPA 



We base our discussion on the Lagrangian of the relativistic mean field (RMF) theory: 

Crmf = 4> N (ijfid' 1 - M N - g a a - g UJ j fl uj t * - g P ^ fi T a p afJ *)ip N 
+ \d ll ad^a - ^mla 2 - ^ g 2 a 3 - ^g 3 a 4 

- \G%G^ + \ m 2 p plp^ , (2.23) 

where the standard notations are used: ipN denotes a baryonic field of mass Mm- er, lo^ and p° are a, lo and p meson 
fields with m a — 511MeV, m w = 783MeV and m p — 770MeV, respectively. F M „ and G" are the antisymmetric field 
tensors for ui and p mesons. The constants g a = 10.0, g u = 12.6 and g p — 4.6 are the coupling constants for the 
interactions between mesons and nucleons. g 2 = — 7.2fm _1 , g 3 = 0.62 and C3 = 71.3 are the self-coupling constants for 
a and w mesons. The quoted numbers are taken from the TM1 parameter set published by Sugahara & Toki [|3l[(see 
also H). 

The approximate self-energy employed in the mean field theory is given as follows. For simplicity, only the con- 
tribution from er meson, that is the scalar part, is shown. The extension to other meson contributions is done in an 
analogous way. 
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E,(l',l) = g„ 5(1-1') ■ a(l) 



(2.24) 



where a meson is assumed to be a classical field and obeys the equation of motion: 

(d v d v + ml) a + g 2 a 2 + g 3 cr 3 = - A^ 1 a + g 2 a 2 + g 3 a 3 = - g^tf^Jf] 



(2.25) 



Here A^" 1 is the inverse of the free propagator. Solving Eqs. ( 2.24 ), ( 2.25 ) with the Dyson equation of nucleon 
consistently for stationary and uniform matter, we obtain RMF for the nuclear EOS. 

as 



On the other hand, using the same self-energy and Eq. ( 2.25 ), we can evaluate S defined by Eq. ( 2.1 



5(1', 1,2', 2) = g 2 a 5(1' -1)5(2' -2) ■[- A- 1 (1-2) +2 g 2 a + 3 g 3 a 2 ] 5(1 - 2) , 



(2.26) 



where a is the expectation value of a in th e sta tic uniform matter. It is obvious that this approximation is conserving 
in the above sense. Putting this into Eq. (2.16), we get the equation for L(l, 1 + , 2, 2 + ) as 



1+2,2+) = L (l,l+,2,2+) - dx 3 dx 4 L (l,l+,3,3 



A-!(3-4) 



L(4,4+,2,2+) 



(2.27) 



Here A = Aq" — 2g 2 a — 3 g 3 W . This is well known RPA except for the modification of the meson propagator 
due to its self-coupling. Thus we established the thermodynamic consistency between RMF and RPA. 

It should be noted that the modified meson propagator is also obtained by expanding the Lagrangian around the 



Since the stationary solution is nothing but the 
ation around the mean field in the above approximation. 



stationary point up to the quadratic order of perturbation 
mean field, the correlation is described by the harmonic oscilL 
In this sense, it is natural that we get RPA consistent with RMF in the static and long wave length limit. 

The above argument is in a rigorous sense applicable only to the conserved quantity such as the baryonic number. 
In general the axial vector part of the nucleon weak current is not a conserved current. However, we extend the above 
method to the axial vector contribution, since we can thus argue both vector and axial vector contributions on the 
same basis. Moreover, the interpretation that we describe the correlations by the harmonic oscillations around the 
static and uniform mean field is still valid. Thus the basic equations of RPA are summarized as 



Tr[T a L(\ ~ 2)T b ] = Tr[T a L (l - 2) T b ] 

fdx 3 dx 4 Y / Tr[T a Lo(l-3)T c } ■ V p c ot (3 - 4) • Tr [T c L(4 - 2) T b ] 



(2.28) 



Here L(l — 2) and Lq(1 — 2) are the abbreviations of L(l,l+,2,2+) and Lo(l, 1 + , 2, 2 + ), respectively. T's denote 
gamma matrices, 1,7^,7^75, and the matrix structure of the response functions is explicitly indicated by the trace 
operations. V po t is the nuclear potential mediated by mesons and is given by the modified propagators as 



V p s ot (l-2) = -^A(l-2) 



(2.29) 



for a meson, for example. The response function L(l — 2) is analytically continued to the retarded Green function, 



which in turn gives the structure functions by Eq. (2.13). It is also possible to work in the real time formalism from 
the beginning, which we did in this paper. Regardless, they are equivalent to each other f48| . The explicit expressions 
of Tr[T a L (l - 3) r c ] are given in Appendix. 

It is clear from the above derivation that RPA is an approximation obtained from the self-energy evaluated up to 
the first order of the coupling constant and that processes such as collisions of nucleons are ignored. This can be 
understood from the fact that RPA is also obtained from the collisionless Boltzmann equation |3^]. The collisional 
effect might be important particularly in the low density regime [ p4]p5| : |2"7jl a nd could in principle be included in the 
present formulation by taking higher order corrections of the self-energy p9fl. It is, however, very difficult to do this 
in practice. In this paper we make more of the consistency of the approximations used for reaction rates and EOS we 
have currently at hand. The study of possible effects of scatterings which was discussed in the introduction will be 
published elsewhere pq |. 

In RMF we ignore the negative energy contribution to the nucleon spinor. In so doing, the p ositiv e energ y part 
is defined in a density-dependent way. This extra density dependence should have made Eqs. ( 2.27 ), ( 2.28 ) more 
complicated, which we ignored in this paper. It should be noted that this inconsistency is quite minor except for 
small parametr regions very close to the phase boundary with which we are not concerned here. 
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C. Residual Interactions 



The Lagrangian shown above does not include the contributions of the pion and the tensor coupling of the p meson 
since they do not contribute to the mean field |30|]3l| l. To describe fluctuations in spin and isospin around the mean 
field these residual interactions should be added to the above Lagrangian as pil 



C — C 



£ P = 



- V 



N 



9tt 

2M N 



7 M 75^7r Q 



C„ 



2M 



N 



(2.30) 



with gl/An = 14.08, ^/4tt = 0.41 and C p = 6.1. Accordingly is added to Fs in Eq. ( |2.28| ). 

It is also known that there are short ranged repulsive correlations for these isovector channels. They are conveniently 
included by the phenomcnological Landau - Migdal parameters |52|. In this paper they are implemented, following 



Horowitz & Piekarewicz []34 57 1 , by modifying the meson propagators as 



9 9 



(I 2 



(2.31) 
(2.32) 



with the Landau - Migdal parameters g'^ = 0.70 and g' p = 0.30. Unlike the conventional non-relativistic treatment, 
we take different parameter value for p mesons from pions. which better reproduces the electron scattering data as 
well as the spin-transfer observables in (p, n) reactions |34|,[35|,^1 -53 57[|. 

As for the electromagnetic interaction, we make the Thomas-Fermi approximation, where the photon propagator is 
replaced by the screened one as 



Itf 



(2.33) 



with q^ F = 4 e 2 7r 1//3 (3 n e ) 2 ^ 3 , w here e and n e are the electron charge and number density, respectively. Then it was 
added to the y pot 's in Eq. fl2.28| ). 

It is finally noted that since the neutrino couples only to the spin-transverse correlations in RPA, the pion (spin- 
longitudinal correlations) does not contribute to the neutrino reaction rates considered here. 



III. RESULTS 



A. Neutral Current Reactions 



In this section, we discuss correlation effects on the neutral current reactions, that is, neutrino - nucleon scatterings. 
The high density regime, p\> > 2 x 10 1 4 g/c m 3 , and the low density regime, pi, < 2 x 10 14 g/cm 3 , are discussed separately, 
since the behaviors of R\{k) in Eq. ( |2.8| ), which is refered to as the vector current correlations in the following, are 
qualitatively different between the two regimes. 

The vector current correlation is mainly induced by the scalar isoscalar meson a and the vector isoscalar meson u, 
while the axial vector current, which is more important for the neutrino opacity, is dominantly affected by the vector 
isovector meson p. 

As shown shortly, the structure functions of nucleons are in general narrow peaked as a function of the transferred 
energy in the sense that their width is much smaller than the transferred three momentum. This feature reflects the 
fact that the nucleon mass is typically much greater than the neutrino energy and the scattering is almost isoenergetic. 
If we assume that the scatterings are exactly isoenergetic, that is, the structure functions are proportional to the delta 
function of the transferred energy, then we obtain the well known Bruenn's formula Q which we use as a reference in 
this paper. In some cases this approximation alone overestimates the cross section significantly |54j . This issue will 
be addressed again in the following sections. 
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1. high density regime 



First we consider the high density regime, p& > 2 x 10 14 g/cm 3 . Fig. [i] shows the structure functions i?i(fco, | A^l) 
and i?2(^0i | fe |) as a function of the transferred energy fco for the density, temperature and proton fraction given in 
the figure. The transferred three momentum | k | is taken to be a typical neutrino energy ~ 3 T. The long dashed 
curves correspond to the structure functions of the non-interacting nucleons. The width of the structure functions 
gives a measure of average energy exchange between neutrino and nucleon. As stated above, we obtain the Bruenn's 
formula [Q if we neglect this width of these functions entirely and approximate the structure functions by the delta 
functions of fco as: 

R 1 {k , \t\)^2wS(k ) ■ [h v 2 r)n + h p v 2 Vp ] (3.1) 
R 2 {k ,\t\)^2wS(k ) ■ [v\ 2 Vn + h p A 2 r) p ] (3.2) 

with tin — / - — r-r /at(p)[1 — /aKp)] f° r N = n,p. Here /at(p) is Fermi-Dirac distribution function for the nucleon. 
J (2tt) j 

Note that the long wave length limit | fc | ^ is also taken in this approximation. It is obvious from the figure that 
this approximation is not very good for this large momentum transfer which is common in the hot proto neutron star 
It should be noted that the neglection of the energy transfer leads to the overestimation of the reaction rate. 



The short dashed curves are obtained from the first term including Lq on the right hand side of Eq. (2.28), where 
only the modified dispersion relation in medium, that is, the effective mass and potential of nucleons is taken into 
account. Note the effective mass of the nucleon becomes as small as ~ 0.6Mn in RMF (see The low effective 

mass renders the structure functions a little bit wider than for the non-interacting case, since the energy exchange 
between neutrino and nucleon is more facilitated, and lowers their amplitudes as well. It is important to note that 
with this modification alone, the neutrino scattering rates are considerably reduced However, as shown in the 
preceding section, the inclusion of the effective mass alone is not consistent with EOS obtained in RMF. 

The solid curves in the figure shows the structure functions with the RPA correlation. As already mentioned, the 
scalar isoscalar meson a and the vector isoscalar meson to are the dominant agents of the correlation for i?i(fc). In the 
high density regime, the contribution of lo dominates over that of a, making the nuclear force effectively repulsive and 
having the structure function even smaller for small fco. O n the other hand, the p meson is the dominant mediator of 
the nuclear correlation for R 2 (k). Since the typical transferred momentum is still small enough for the short ranged 
repulsive force described by the Landau - Migdal parameters to be dominant, the response of this channel is also 
suppressed. The feature mentioned here is common to the structure functions in this density regime. 

As the temperature decreases, the structure function is more skewed toward the positive energy transfer fco > 0, 
reflecting the fact that the extraction of energy from the nuclear medium becomes more difficult due to the Fermi 
blocking for the down scattered nucleon. This is shown in Fig. |2|. 

On the other hand, as the density becomes higher, the nucleon effective mass gets smaller (M^ ~ 0.4Mjv at 
pt = 5 x 10 14 g/cm 3 ), making the amplitudes of the structure functions, both for R\(k) and R 2 (k) 1 smaller without 
RPA. The oj meson becomes even more dominant over the a meson and thus RPA further reduces the amplitude of 
i?i(fc), as shown in Fig. ||. 

The case of the proton fraction Y p — 0.1 is shown in Fig. |]. The effect of RPA increases with decreasing Y p , 
which could be understood as follows. The vector current correlation comes mainly from the neutron sector, since the 
coupling constant for proton hy ~ 0. Since the neutron density increases as the proton fraction decreases, the vector 
current correlation becomes greater, reducing the structure function R\(k). 

The total scattering rate of the neutrino with the incident energy E™ is given by the integration of the structure 
function with respect to the transferred momentum k and energy fco: 
■ d 3 q out 1 



p Ji put 1 

R t0t (K n ) =J-^s 2E ou t R( E ™ K u \ cosfl) [1 -/„(£?*)] 
1 f°° n , „ /- fc o mQX E'l n -k 1 



(27T) 



poo rK Q jpzn h. i 

Jo L k \e?2E?* R{E " 1 ' E ° Ut > C0S9) [1 _ MKUt)] ' (3 - 3) 



with fco™" 21 = mi n(k , 2E™ — k ). If we ignore the energy exchange and the Fermi blocking for the scattered neutrino, 
and insert Eqs. (3.1) and ( |3.2| ) into Ri and R 2 in the above equation, we obtain the standard approximate formula Q 



Rtot[ET) w £ G% (Kl)2 |^ j + 3 (/4)2 j | (3 4) 

N 
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which is frequently used in literatures. In the following we use this rate as a reference in evaluating quantitatively the 
suppression factors of the total scattering rates due to the correlation. It should be noted, however, that the difference 
between these rates also comes from the inclusion of the effective nucleon mass and the finite energy and momentum 
transfer, which are neglected in Eq. (3.4). 

Fig. U shows the ratio of R tot evaluated for the first term in Eq. ( [2.2S| ) by using the exact formula given in Eq. (3.3) 
to R tot by using the approximate expression given in Eq. (3.4) as a contour map in the density and temperature plane 
for the given proton fraction. The contributions from R\ and i?2 are shown separately. The incident neutrino energy 
is assumed to be 3 T, and the blocking factor was dropped again. 

It is again clear that the reaction rates are substantially smaller than the standard ones before RPA is included. 
This is mainly due to the decrease of the nucleon effective mass. In fact, as the density increases, the reduction 
becomes greater, which is common to the vector current part from Ri and the axial vector current part from R 2 . 
The non-monotonic temperature dependence is understood as follows: In the standard approximation formula, as 
mentioned above, the energy exchange between neutrino and nucleon is neglected. In the low temperature regime, 
this neglection tends to overestimate the Fermi blocking for the scattered nucleon, reducing the standard rate. As 
a result, the suppression factor becomes larger as the temperature increases. For even higher temperature, however, 
the width of the structure function cannot be ignored and lowers the total scattering rate just like the neutrino - 
electron scattering. Since the latter effect surpasses the former at some temperature, the suppression factor begins to 
decrease again above this temperature, which becomes larger as the density increases. These features were discussed 
in the paper by Schinder H] where he studied the effect of the recoil of nucleons. 

Fig. ^ shows the suppression factor for the reaction rates with the inclusion of RPA. It is evident that both the 
vector current and axial vector current contributions are further suppressed, which is expected from the above results 
for the structure functions. Indeed the reaction rates become less than half the standard rates around the saturation 
density. Since the delayed explosion is very sensitive to the neutrino luminosity and the energy, this suppression could 
have an significant influence on the final outcome of the core collapse. 

we can understand the different density and temperature dependences of the vector and axial vector contributions. 
In the former case, the RPA suppression is strongly density dependent, since the nuclear force is determined by the 
competition of the attractive a meson and the repulsive to meson. As a result, the reaction rates are more strongly 
suppressed in the higher density regime. On the other hand, the balance between the attractive p meson and the 
repulsive short range force described by the Landau - Migdal parameter is determined by the momentum transfer, 
thus more temperature dependent. Since the nuclear force is more repulsive for the smaller momentum transfer, that 
is, in the low temperature regime, the rate is reduced further there. 

Figs. and || shows the suppression factors for the different proton fractions. As expected, the axial vector part is 
more affected, since the isovector p meson is convening its correlation. It is also clear that the suppression is stronger 
for the smaller proton fraction. 



2. low density regime 

The low density regime considered here is characterized by the appearance of nuclei for low temperatures |2^,^9|,^5| . 
Although we cannot treat the coexistence phase of nuclei and nucleons by the current formulation and consider here 
only the gas phase where no nuclei exist, the precursor of the phase transition is already imprinted in this regime. 

On the contrary to the previous section, the a meson dominates the u> meson and the isoscalar nuclear force becomes 
attractive. As a result, the correlation enhances the reaction rate coming from the vector current part. This can be 
seen in Fig. |9|, where the structure functions R\ and Ri are shown for the temperature near the phase boundary. 
The enhancement of R\ in the vicinity of fco = suggests the existence of the static unstable mode [M. In fact, 

Ri diverges on the so-called spinodal line From the thermodynamic consistency discussed above, (j^r) a l so 

diverges on this line. This behavior is better indicated in Fig. [UJ in which the enhancement factor with respect to the 
standard Bruenn's formula is plotted as a contour in the density temperature plane. As the temperature decreases, 
the reaction rate becomes larger. The lower boundary of the computation region roughly corresponds to the spinodal 
line obtained by the current model. Incidentally, the axial vector contribution is hardly affected since the a and u> 
mesons are minor contributors for this channel. 

Since in reality this liquid-gas phase transition occurs not by the spinodal decomposition but by the nucleation, the 
actual phase boundary does not correspond to the spinodal line except for the critical point, and the spinodal line is 
never reached |^,^9|,^5|. The critical temperature obtained from RPA is about 15MeV, which is actually consistent 
with the value obtained from EOS by Shen et al. p8| , p9] |. 

Although the enhancement factor is still relatively large in the region not very close to the phase boundary, this 
effect has only minor importance for the neutrino opacity, because the vector current contribution is typically much 
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smaller (~ 1/4) than the axial vector contribution. Provided that the axial vector part is slightly suppressed in 
this regime due to the repulsive nuclear interaction for that channel, the total scattering rate will not be increased 
significantly unless the critical point is approached very closely. 



B. Charged Current Reactions 



The RPA correlation is calculated for the charged current reactions just the same way as for the neutral current 
reactions except that the vector isovector p meson is the only mediator of the correlation. It is noted again that 
since only the spin-transverse correlation matters for the neutrino reaction rates in RPA, the pion does not make any 
contribution. As a result, for relatively small momentum transfer of our current interest the nuclear force is repulsive 
and the correlation always tends to suppress both the vector current and the axial vector current contribution to the 
charged current reaction rates like the axial vector part of the neutral current. This is true not only for the high 
density regime but also for the low density regime unlike the vector current part of the neutral current reaction rates 
which is enhanced by the liquid - gas phase transition. 

Fig. |ll| shows a typical modification of the structure functions around the saturation density for the symmetric 
nuclear matter. Depicted are Ri(ko, \k\) and i?2(fco, 1 k |) for the e~ + p — > v e + n reaction with the outgoing neutrino 
energy of ~ fi e , the electron chemical potential. It is clear again that the decrease of the nucleon effective mass 
alone (the short dashed curve) accounts for more than half of the suppression with respect to the non-interacting 
case (the long dashed curve). Moreover, as the electron chemical potential is quite large, ~ 250MeV, the average 
energy transfer from the incident electron to the nucleon becomes also large. Hence the complete neglection of energy 
transfer assumed in the approximate formula 0] leads to sizable overestimation of the reaction rates for this high 
energy neutrino emission. 

The results change qualitatively for the asymmetric matter, since the potential for the neutron is different from 
that for the proton in medium. This is understood as follows. In the non-relativistic limit the dispersion relation for 
the nucleon is written as 



k\ 
2M* N 

and the potential difference is given in RMF by 



E N(k) = 1_J_ + up , (3.5) 



MJ^ = Ul ot - Uf* = 2g p (p) = 2^(n p -n n ) . (3.6) 

nip 

Thus, in neutron rich matter, the proton is more strongly bound than the neutron. This difference of potentials 
serves as a threshold for the p — > n reaction as the mass difference of nucleons does for the reaction in vacuum. It is, 
however, ignored in the approximate formula ||, which is given by 

R^ko, \Te\) « 2tt 5(k Q - A) • g v 2 ^ (3.7) 
R 2 (k , » 2ir5(k Q - A) ■ g A 2 Vpn (3.8) 

[2 cfip 

with rj pn — I - — — f p (p)[l — f n (p)]- Hence the structure functions in this approximation are delta functions placed 
J (27r) J 

at fc = A = M n - M p . 

With the potential difference taken into account, however, the structure functions are shifted by —AUg"*, thus 
they do not agree with the approximate ones even in the static and long wave length limit with —> Mjv- This is 
seen in Fig. 12, where the structure functions for the asymmetric matter (Y p = 0.3) are shown for the non-interacting 
case (the long dashed curve) as well as the case with the nucleon effective mass and potential alone included (the 
short dashed curve) and the case with RPA also included (the solid curve) . In the non-interacting case, the structure 
functions are located around fco ~ A, while in the other cases they are moved toward the positive energy transfer as 
stated above. 

The total v e emission rate is given for an outgoing neutrino energy E° ut as 

R totr E out, = / R(El n , EZ U \ COSO) f e (ET) 

J (27r) 3 2E m 2E ' Ve7 v ' / Je\ e J 

1 r°° f^o zpout _[_ h i 

= —— 2irkdk dk v - t ° : R(El n , E° ut , cos 8) fJEl n ) , (3.9) 

(2tt) 3 J J k - ° E° ut 2E™2E° ut V e ' " ' ' JeK e ' ' y ' 
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where k^ = y/m* + (E° ut ± k) 2 — E° ut . Inserting Eqs. (3.7), (3.8) into the above equation, we obtain the standard 
approximate formula: 



R tot {E ou t) „ ^ {gv 2 + 3gA 2 } ^ {E out + A) 2 y 1 _ ( go r+ A ) /-(^ + A ) • ( 3 - 10 ) 

For the symmetric nuclear matter, t he d ensity a nd te mperature dependence of the suppression factor of the total 
emission rates, that is the ratio of Eq. (3.9) to Eq. ( 3.10| ), is essentially the same as th at o f the axial vector contribution 



to the neutral current reaction. For the asymmetric matter, as understood from Eq. ( |3.9| ), the k integration has greater 
contribution from larger momentum transfer due to the shift of the structure functions, which tends to get the reaction 
rates larger thanks to the larger phase volume. However, as E° ut comes closer to /Lt e , the electron population starts 



to deplete f e (Ej n ) < 1 and, as a result the emission rate becomes smaller. This occurs for smaller E° ut in Eq. (3.9) 
than in Eq. ( 3.1 0| ) due to the large width of t he structure function. It follows that the total emission rate is greater 



than the approximate one given by Eq. (3.1C) for neutrinos with energy E° <C fi e and smaller for neutrinos with 
E° ut < p e . This is indicated in Fig. [l|. Thus lower energy electron neutrinos are produced more efficiently than 
predicted from the standard approximate formula. It is also seen that the RPA correction is not very large since the 
nuclear force becomes less repulsive for the large momentum transfer considered here. 

Fig. [l4]displays the structure functions for a lower temperature. It is evident that their widths became smaller. This 
is because nucleons become more degenerate as the temperature is lowered and the large energy transfer is required 
to overcome the difference of chemical potentials between neutron and proton. The comparison of Fig. |l5|, which 
represents the emission rates for the same temperature, with Fig. [l3| shows that the emission of high energy neutrinos 
is severely suppressed for the low temperature, because it is contributed from the high energy tail of degenerate 
electrons. It is also seen from the figure that the approximate formula tends to overestimate the Fermi blocking factor 
for the outgoing neutron by neglecting the finite energy transfer, thus underestimating the emission rate. The RPA 
correction is even smaller in this case, since the narrow structure functions favor higher momentum transfer. 

As the density is increased or the proton fraction is decreased, the potential difference becomes larger. Hence the 
features mentioned above become more prominent, which is demonstrated in Figs. [l6|and In the latter case, since 
the electron chemical potential is also decreased, it is possible for small Y p and T that the structure functions and 



the electron distribution function in the integrand of Eq. (3.9) do not give a significant overlap and the emissivity is 



strongly reduced. This is well known suppression of the URCA process |38| . 

The structure functions for the n — > p reactions are derived from those for the p — > n reactions by using the detailed 
balance relation expressed as 

R {n ' p) {k) = e ^o-^pn) #(p,»)(_fc) , (3.11) 

where A^i pn = fi p — fj, n is the difference of the nucleon chemical potentials. In this paper, however, the n — > p 
structure functions were calculated in the same way as the p — > n counterparts, and we confirmed directly that the 
above relation is actually satisfied. 

As expected, the structure functions for the n->p reactions are shifted toward the negative energy transfer due to 
the potential difference for proton and neutron in asymmetric matter, which is shown in Fig. The RPA correction 
accounts for about half of the total suppression in this case. The total emissivity of T> e via e + + n — > V e +p is obtained 
from Eq. ( |3.9| ) with the following exchanges: p <-> n, e «-» e + and v e «-> v e . The important difference from the v e case 
comes from the fact that positrons are not degenerate in the supernova core. The shift of the structure functions in 
the direction of the negative energy transfer gives greater weights to positrons with lower energies which are more 
populated, leading to larger emissivity of v e . On the other hand, the RPA correction tends to reduce the rate as in 
the v e case. These trends are shown in Fig. [l9], where the V e emissivity is plotted as a function of the energy of the 
emitted neutrino. The potential difference for neutron and proton provides the emitted neutrino with certain energy 
(~ 30MeV in this case). Although the correlation reduces the emissivity substantially, the emissivity of high energy 
neutrinos are larger even after RPA is included than for the approximate formula. 

Thus far, we discussed the v emissivities. The corresponding absorption rates of v e via v e + n — > e + p and 77 e via 
V e + p — > e + + n are obtained by the detailed balance equation: 

R%(E V ) = e /^-M Rl%(E u ), (3.12) 

where the neutrino che mical potential is defined as /i„ e = /i e + ji p — fi n for v e and fj^ r — —fJ-u e f° r v e . This relation is 



easily proven from Eq. (3.11) 
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IV. CONCLUSION 



In view of the application to supernova simulations |59j6fJ, we have calculated the neutrino - nucleon reaction 
rates in hot and dense supernova cores with the RPA correlations taken into account. The approximations are based 
on the Lagrangian used in RMF for the nuclear EOS, where the nonlinear a and to meson terms are included. This 
ensures the thermodynamic consistency in the static and long wave length limit between the EOS we use in the 
supernova simulations and the vector current correlation of nucleons which gives partial contributions to the neutrino 
- nucleon scattering rates. The same method is extended to the axial vector current as well as to the isovector current 
correlations with the residual interactions added which are described by the tensor coupling of the p meson and the 
so-called Landau-Migdal parameters. 

We have found that the neutral current reactions are suppressed substantially in the high density regime, p > 
10 14 g/cm 3 , due to the repulsive nature of the nuclear forces. In the lower density regime, however, the vector current 
contributions are enhanced as the temperature is decreased and the boundary of the gas phase, which is placed by RPA 
consistently with EOS, is approached. Although this is interesting itself, the total scattering rates are not enhanced 
significantly except for the very vicinity of the critical point, since the axial vector contribution is dominant over the 
vector current one and is unaffected by the liquid - gas phase transition. We have also shown that the neglection of 
the finite energy exchange between neutrino and nucleon, which is assumed in the standard approximation formula, 
can lead to a significant overestimation of the total reaction rates in the high temperature regime where the typical 
incident neutrino energy is no longer much smaller than the nucleon effective mass. 

In the case of the charged current reactions, it is more important to take proper account of the potential difference 
for neutron and proton than to consider the RPA correlation, since it serves as a threshold of the reaction just as a 
mass difference does in vacuum. In fact, with the potential difference included the v e emissivity via e + p — > v e + n is 
much smaller for higher energy neutrinos than expected from the standard formula. On the other hand, the emission 
of low energy antineutrinos via e + + n — > v e + p is reduced. In both cases, the RPA correlation reduces the total rates 
further. 

These modifications of the neutrino nucleon reaction rates could have an important impact on the mechanism 
of collapse-driven supernovae, since they are the dominant opacity sources for neutrinos. The smaller opacity leads 
to the larger neutrino luminosity, hence to the greater heating rate behind the stagnated shock wave, which gives a 
better odds to get a successful explosion. The changes of neutrino mean free paths, on the other hand, could also 
affect the strength and region of convections in the proto neutron star, which is one of the key ingredients to determine 
the neutrino luminosity and energy. These issues will be discussed in detail elsewhere J6l| . 

It is clear that RPA does not include everything. As understood from the fact that RPA is consistent with the 
mean field theory for the equation of state, RPA studies the response of the mean field to the external disturbance. 
The additional correlations might be induced by some reactions such as collisions of two nucleons via spin and isospin 
dependent reactions, as pointed out by some authors. Although these effects can in principle be included in the above 
formulation by taking the corresponding higher order corrections to the self-energy, it is extremely difficult to achieve 
it in practice. In this paper we employed RPA respecting the consistency between the calculated reaction rates and 
EOS we have currently at our disposal. Nevertheless, these issues are remaining to be studied further. 
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APPENDIX: CAUSAL GREEN FUNCTIONS 



Once the formulation is established in the imaginary time formalism, it is more convenient to wor k wit h quantities 
which are analytically continued from the imaginary time domian to the real time domain. Then Eq. ( 2.2S| ) is still valid 
with the response functions replaced with the retarded Green functions. These rea l tim e retarded Green fun ction s 
Eq. ( 2.12| ) are derived from the corresponding causal Green functi ons d efined by Eq. ( 2.11 ) via the relation Eq. ( 2.13 ). 
Thus all we have to do is to calculate II™ and then to solve Eq. ( p. 28 ) which is now rewritten in momentum space as 



'(*) • Koti.k) 



(Al) 
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Here II™ ab (fc) is derived by Eq. (2.13) from n^^fe) which is expressed as 



ill% ab (k) 



Tr 



d 4 q 
(2^4 



T a G{q + k)T b G(q) 



where T's stand for 1, 7^, 7^75. G(fc) is a single particle Green function of nucleon given by RMF: 
G{k) = G F (k) + G D (k) 



G F (k) 



k 2 - M 



N 



G D (k) = $ + M* N ) ■ 2mS(k 2 - M* N 2 ) {Q(k ) f N (k ) + 6(-fe ) 



(A2) 

(A3) 
(A4) 

(A5) 



where k = (k — Sy, fe), Af^ = Af^r + Eg with Eg, Ey the sclar and vector nucleon self-energies given by RMF, and 
In , /jv are Fermi distribution functions for nucleon and antinucleon, respectively. Note that since the above Green 
function G(k) already includes the effective mass and potential, n^(fc) is different from the response function of 
non-interacting nucleons. As shown in the above equation, the Green function is decomposed into two contributions, 
that is, the Feynman part Gp(k) and the density-dependent part Gi)(fc). As a usual practice, we ignore the vacuum 
polarization coming from the Feynman part alone in calculating I\° N {k). Taking the proper linear combination of 



ab 



we obtain the retarded Green function for the nucleon weak current, from whic h u sing again Eq. (2.13), we 
finally arrive at the structure functions of nucleons. The explicit expressions of Eq. (A2) are found, for example, 
in I2IJ2I& 
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FIG. 1. The vector and axial vector structure functions Ri(ko, | fc |) and i?2(fco, | k |) for the neutral current as functions 
of the transfer energy ko. The density pj,, temperature T, proton fraction Y p and the absolute value of the transferred three 
momentum | k | are shown in the figure. The long dashed curve represents the non-interacting case, while the short dashed 
curve shows the case for the effective mass of the nucleon being taken into account. The solid curve shows the results of RPA. 



FIG. 2. The vector and axial vector structure functions R\(ko, | k |) and i?2(feo, | fc |) for the neutral current as functions of 
the transfer energy ko. The density pt and proton fraction Y p are the same as in Fig. |l[ but the temperature T is smaller in 
this case, and the transferred three momentum | fe [ is scaled as ~ 3T. The explicit values are shown in the figure. 



FIG. 3. The vector and axial vector structure functions i?i(feq> | k ) and i?2(fco, | fe |) for the neutral current as functions of 
the transfer energy fco. The parameters are the same as in Fig. hj but for the higher density, pt, = 5 X 10 14 g/cm 3 as indicated 
in the figure. 



FIG. 4. The vector and axial vector structure functions RUko, \ k |) and R2(ko, \ k |) for the neutral current as functions of 
the transfer energy ko. The parameters are the same as Fig. hi but for the different proton fraction Y p = 0.1. 
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FIG. 5. The suppression factors for the total scattering rates, R tot {E]J l ) defined in Eq. (1-3), from the expressions of Bruenn 
for vector currents (upper figure) and for axial vector currents (lower figure). Here E™ is taken as El" — 3T. The proton 
fraction is chosen as Y p — 0.3. The suppressions are caused by the effective mass for nucleons and the effect of the transfer 
energy spectra instead of the 5-function. 



FIG. 6. The suppression factors for the total scattering rates, R {E^J 1 ) defined in Eq. (3.2), from the expressions of Bruenn 
for vector currents (upper figure) and for axial vector currents (lower figure). E™ = 3T and Y p = 0.3 are used here also as in 
Fig. p. In this case, the suppressions are caused by the RPA correlations in addition to the effects considered in Fig. |E| 



FIG. 7. The suppression factors for the total scattering rates, R tot (El n ) defined in Eq. (3.2), from the expressions of Bruenn 
for vector currents (upper figure) and for axial vector currents (lower figure). The conditions are the same as those in Fig. [], 
but the proton fraction is reduced here to Y p — 0.1. 



FIG. 8. The suppression factors for the total scattering rates, R tot {E l J l ) defined in Eq. (3.2), from the expressions of Bruenn 
for vector currents (upper figure) and for axial vector currents (lower figure). The conditions are the same as those in Fig. [], 
but the proton fraction is changed here to Y p = 0.5. 



FIG. 9. The structure functions R\(ko, | fc |) and i?2(fco, I k I) f° r the neutral current in the low density regime. The density 
pb, temperature T, proton fraction Y p and the absolute value of the transferred three momentum | fc \ are shown in the figure. 
Only the effective mass of the nucleon is taken into account for the short dashed curve, while RPA is included for the solid 
curve. 



FIG. 10. The enhancement factor for the vector current contribution (upper) and the suppression factor for the axial vector 
current part (lower) for R tot (EU 1 ') defined in Eq. ( 3.3) are shown for the low density regime. Here the proton fraction is Y p — 0.3. 
The dark regions are outside the computational grid and roughly correspond to the coexistence phase of nucleons and nuclei, 
where a special care is required to treat the non-uniform matter. 



FIG. 11. The structure functions Ri(ko, | fc |) and i?2(fco, [ k \) for the charged current reactions, where protons change into 
neutrons. The density pi,, temperature T, proton fraction Y p and the absolute value of the transferred three momentum | k 
are shown in the figure. The long dashed curve represents the non-interacting case while the short dashed curve shows the case 
for the effective mass of the nucleon being taken into account. The solid curve shows the results of RPA. 



FIG. 12. The structure functions R\(ko, \ k |) and R2{ko, \ k |) for the charged current reactions, where protons change into 
neutrons, for the proton fraction, Y p = 0.3. The three momentum transfer | fe | is also changed to the electron chemical potential. 
The notations for curves are the same as in Fig. |ll]. 



FIG. 13. The total i/ e emission rate defined in Eq. (3.9) as a function of the emitted neutrino energy E° ut . The upper panel 
shows the rate coming from the vector current and the lower panel represents the axial vector current contribution. The long 
dashed curves are obtained from the Bruenn's approximation formula, while the short dashed curves represent the results with 
the nucleon effective mass and the isovector potential being included. The solid curves correpond to the RPA results. 



FIG. 14. The structure functions Ri(ko, \ k |) and Rz(ko, \ k |) for the charged current reactions, where protons change into 
neutrons. The density pt, temperature T, proton fraction Y p and the absolute value of the transferred three momentum | fe 
are shown in the figure. The long dashed curve represents the non-interacting case while the short dashed curve shows the case 
for the effective mass of the nucleon being taken into account. The solid curve shows the results of RPA. Here the temperature 
is T = lOMeV. The three momentum transfer | fc | is also changed to the electron chemical potential. The notations for curves 
are the same as in Fig. [ll]. 
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FIG. 15. The total emission rate denned in Eq. (3.9) as a function of the emitted neutrino energy E° ut . The upper panel 
shows the rate coming from the vector current and the lower panel represents the axial vector current contribution. The long 
dashed curves are obtained from the Bruenn's approximation formula, while the short dashed curves represent the results with 
the nucleon effective mass and the isovector potential being included. The solid curves correpond to the RPA results. Here the 
temperature is T = lOMeV. 



FIG. 16. The total u e emission rate defined in Eq. (3.9) as a function of the emitted neutrino energy E° ut . The upper panel 
shows the rate coming from the vector current and the lower panel represents the axial vector current contribution. The long 
dashed curves are obtained from the Bruenn's approximation formula, while the short dashed curves represent the results with 
the nucleon effective mass and the isovector potential being included. The solid curves correpond to the RPA results. Here the 
density is pt — 5 x 10 14 g/cm 3 



FIG. 17. The total z/ e emission rate defined in Eq. (3.9) as a function of the emitted neutrino energy E° ut . The upper panel 
shows the rate coming from the vector current and the lower panel represents the axial vector current contribution. The long 
dashed curves are obtained from the Bruenn's approximation formula, while the short dashed curves represent the results with 
the nucleon effective mass and the isovector potential being included. The solid curves correpond to the RPA results. Here the 
proton fraction is Y p — 0.1. 



FIG. 18. The structure functions Ri(ko, | k |) and J?2(fco, | k |) for the charged current reactions, where neutrons change into 
protons. The density pi,, temperature T, proton fraction Y p and the absolute value of the transferred three momentum |fc| arc 
shown in the figure. The long dashed curve represents the non-interacting case while the short dashed curve shows the case 
for the effective mass and the isovector potential of the nucleon being taken into account. The solid curve shows the results of 
RPA. 



FIG. 19. The total 77 e emission rate as a function of the emitted neutrino energy E° ut . The upper panel shows the rate 
coming from the vector current and the lower panel represents the axial vector current contribution. The long dashed curves 
are obtained from the Bruenn's approximation formula, the short dashed curves represent the results with the nucleon effective 
mass and the isovector potential being included. The solid curves display the RPA results. 
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